Softening of the insulating phase near Tc for the photo-induced insulator-to-metal 

phase transition in vanadium dioxide 



O 

o 
o 

o 



D. J. Hilton,! B 

R. P. Prasankumar,! S. Fourmaux,'^ A. Cavallcri.'^ D. Brassard,^ 
M. A. El Khakani,2 J. C. Kieffer,^ A. J. Taylor,^ and R. D. Avcritti'0 

^Center for Integrated Nanotechnologies, MS Kill, 
Los Alamos National Laboratory, Los Alamos, NM 87545, USA 
^Universite du Quebec, INRS-Energie et Materiaux et Telecommunications, Varennes, Quebec, Canada, J3X 1S2 
'"^ Department of Physics, Clarendon Laboratory, University of Oxford, Parks Rd. Oxford, 0X1 3PU, United Kingdom 

(Dated: February 1, 2008) 

We use optical-pump terahertz-probe spectroscopy to investigate the near-threshold behavior of 
the photoinduced insulator-to-metal (IM) transition in vanadium dioxide thin films. Upon approach- 
ing Tc a reduction in the fluence required to drive the IM transition is observed, consistent with a 
softening of the insulating state due to an increasing metallic volume fraction (below the percolation 
limit). This phase coexistence facilitates the growth of a homogeneous metallic conducting phase 
following superheating via photoexcitation. A simple dynamic model using Bruggeman effective 
medium theory describes the observed initial condition sensitivity. 

PACS numbers: (71.30.+h), (78.47.-l-p), and (72.80.Ga) 
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Strongly correlated electron materials are character- 
ized by extreme sensitivity to external stimuli which re- 
sult from the subtle interplay between many degrees of 
freedom with comparable energy scales (lattice, orbital, 
electronic and magnetic). High temperature supercon- 
ductivity, colossal magnetoresistance, ferromagnetism, 
and insulator-to-metal transitions are some of the best 
known examples of the exotic behavior exhibited by these 
materials. 

During the past several years, time-resolved optical 
studies have been utilized to investigate dynamics asso- 
ciated with the competing degrees of freedom in these 
materials [H . In many of these studies the goal is to in- 
terrogate the dynamics within a particular phase. How- 
ever, photoinduced phase transitions provide an impor- 
tant complementary approach to investigate the physical 
pathway connecting different correlated electron states as 
well as their mutual competition Q]. 

A photoinduced phase transition can arise following 
impulsive heating or photo-doping and provides a means 
of controUingthe overall phase of a solid on an ultrafast 
timescale . One model system used to study insu- 

lator to metal transitions in correlated electron systems 
is vanadium dioxide (VO2). This material undergoes an 
insulator to metal transition when heated above 340 K 
accompanied by a structural distortion from a monoclinic 
to a rutile phased, S 9- Time-resolved optical studies 
can provide insight into its origin and its technological 
potential[i,[l3,[ilj. 

The interplay between band and Mott-insulating be- 
havior in VO2 has been addressed in the time domain, 
where a limiting structural timescale for the creation of 
the metallic phase suggests the importance of band in- 
sulating character for VO2 Advances in ultrafast 
technology have also enabled more comprehensive inves- 
tigations involving, most recently, femtosecond x-ray ab- 
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sorption spectroscopy |13l . 

In this Letter, we measure the time-dependent con- 
ductivity of VO2 during a photoinduced insulator-metal 
transition. Our dynamic experiments focus on the poorly 
understood near-threshold behavior, where phase sepa- 
ration and domain formation inhibit the measurement 
of conventional transport properties since a macroscopic 
conductivity pathway is not established As a func- 
tion of increasing initial temperature (up to ~ Tc) we 
observe a reduction in the deposited energy required to 
drive the IM transition. Such a response indicates an 
initial condition sensitivity which we interpret as a soft- 
ening of the insulating state due to the existence of metal- 
lic precursors. These precursors facilitate the growth of 
a homogeneous metallic conducting phase following su- 
perheating. A simple dynamic model using Bruggeman 
effective medium theory describes the observed response. 
We emphasize that this percolation model is valid near 
threshold whereas, at hi gher fluences, the transition is 
prompt and nonthermal [l2l . [isl . 



14| . Our results may 



be relevant in other systems which exhibit tendencies to- 
wards phase separation such as spin-crossover complexes 



Our (Oil) oriented VO2 thin film was grown on a (100) 
MgO substrate by radio-frequency (13.56 MHz) mag- 
netron sputtering of a vanadium target (99.99% purity) 
in a Ar and O2 mixture at a pressure of 2 mTorr (cham- 
ber basepressure of 2 x 10~^ Torr) and a temperature of 
500°C [17|. The stability of the deposition process and 
the thickness of the films (100 nm) were monitored in-situ 
by a calibrated microbalance and independently verified 
using standard ellipsometric techniques. The crystalline 
quality, grain size, and orientation of the VO2 films were 
characterized by x-ray diffraction. 

First, we determined the time-independent THz con- 
ductivity of our VO2 films as a function of sample tem- 
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perature. We use the output of a 50 fs titanium:sapphire 
amplifier to generate nearly-single-cycle THz pulses via 
optical rectification in ZnTe and employ standard ter- 
ahertz time-domain spectrometry to measure the THz 
frequency conductivity; further details can be found in 
refs.lisl andflil. In Figure [TJa), we show the transmitted 
THz waveform at 340 K, where the top waveform (dis- 
placed for clarity) results when the sample is heated from 
300 K to 340 K and the bottom results when the sample 
is cooled from 370 K to 340 K. From the time domain 
data, the calculated conductivity of the VO2 film from 
0.25 THz to 2.25 THz is determined (Figure lllb)). 

Figure [Ijc) shows the magnitude of the real conduc- 
tivity as a function of temperature. We also plot the 
temperature-dependent DC electrical conductance mea- 
sured (on the same film) between two metal contacts, 
which shows good agreement with the temperature de- 
pendence of the THz conductivity. The hysteresis is con- 
sistent with the first order phase transition in VO2 and 
the narrow width (< 5 K) and magnitude of the metal- 
lic state conductivity attest to the quality of our films 

In order to study the conductivity dynamics of the 
photoinduccd phase transition, we excite the sample us- 
ing 1.55 cV pulses from the amplified laser and monitor 
the change in the THz transmission as a function of the 
relative delay between the pump pulse and the THz probe 
pulse This film is approximately one optical absorp- 
tion length thick at 800 nm, which results in a nonuni- 
form excitation profile along the direction of propagation; 
however, this does not significantly influence our results. 
We have verified that the induced change in conductiv- 
ity is frequency independent across the bandwidth of our 
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FIG. 1: (color online) (a) THz waveform transmitted tlirougli 
the sample when fieated to 340 K from 300 K (T^ ) and when 
cooled from 370 K (Tj), offset for clarity (b) The optical 
conductivity vs frequency as determined from the transmit- 
ted waveforms in (a) and at 300 K. (c) The temperature- 
dependent hysteresis of the THz conductivity (left scale, 
♦ for T^, • for T^) and of the conductance measured using 

standard DC electrical techniques (right scale, for Tf, 

- for Ti). 
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FIG. 2: (color online) (a) Photoinduced conductivity change 
at 300 K for various fluences where o"(27) is the conductivity 
in the full metallic state (b) Magnitude of the conductivity 
change as a function of fluence at 300 K. 



THz pulse indicating uniform excitation across the THz 
probe beam. 

Figure [IKa) shows the time-resolved conductivity as 
a function of pump fluence at 300 K. There is a rise 
time of '--^lOOps to obtain the conductivity of the product 
metallic phase at all fluences. This is significantly longer 
than the initial excitation pulse, meaning that photocon- 
ductivity arising from carriers excited to the conduction 
band (e.g as occurs in GaAs) is not responsible for the in- 
duced response. Additionally, lattice heating via carrier 
thermalization occurs in approximately one picosecond 
indicating that the conductivity dynamics arc more com- 
plex than simple heating above Tc- Crucially, however, 
for fluences greater than ^^10 mJ cm~^, the deposited en- 
ergy density ('^ 500 J cm"-^) is considerably above what 
is required to heat above Tc (~ 200 J cm~'^ at 300 K). 

Figure [2)Jb) shows the maximum induced conductivity, 
Affmax J at each pump fluence extracted from the data in 
Figure[2Ka). There is a marked decrease in the maximum 
obtainable conductivity with decreasing fluence. At a flu- 
ence of 19.2 mJ cm~^, the induced conductivity at long 
times is 100% of the metallic phase conductivity (i.e. at 
370 K). Extrapolation of the photoinduced conductivity 
change at 300 K, as shown in Fig. 2(b), yields a non- 
zero fluence threshold of ~7 mJ cm~^. The existence 
of a fluence threshold is a well-known feature in photo- 
induced phase transitions, where the cooperative nature 
of the dynamics results in a strongly non-linear conver- 
sion efficiency as a function of the number of absorbed 
photons. In the present case, photoexcitation leads to a 
rapid increase of the lattice temperature, initiating the 
nucleation and growth of metallic domains which coa- 
lesce (i.e. percolate) to yield a macroscopic conductivity 
response. 
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FIG. 4: (color online) The volume fraction, fm, responsible 
for the observed conductivity shown in Figure [!](♦ for Tf and 
• for Ti) and the extracted volume fraction in the insulating 
state (t) as determined from the dynamic BEMT described 
in the text. 



FIG. 3: (color online) (a) Induced conductivity change as a 
function of time at a fluence of 12.8 mj cm~'^ for various ini- 
tial temperatures, (b) Magnitude of the induced conductivity 
(■) and the maximum possible conductivity change (•). The 
black lines are a fit as described in the text, (c) Fluence 
threshold as a function of base temperature. 



We have also measured the photoinduced terahertz 
conductivity as a function of base temperature. An opti- 
cal pump fluence of 12.8 mJ cm~^ was used, less than is 
required to drive the full metallic transition at room tem- 
perature (though, as described above, more than enough 
to superheat the film above Tc). The conductivity dy- 
namics as a function of time at initial temperatures of 
300, 310, and 320 K are displayed in Figure [3ja). From 
Figure E^a) and (b) it is evident that, below Tc, the 
photoinduced change in the conductivity is less than 
the maximum possible induced change. However, with 
increasing initial temperature the induced conductivity 
change increases and obtains the maximum possible value 
at 320 K. Thus, in the insulating phase, there is a de- 
crease in the threshold to drive the sample metallic with 
increasing temperature. At temperatures greater than 
330 K, the photoinduced change in conductivity follows 
the maximum possible induced change (Figure EJb)). 
This occurs since the incident fluence of 12.8 mJ cm^^ 
is sufficient to drive the conductivity to its full metallic 
state consistent with a decreasing threshold. Finally, in 
Figure ^c) we plot the fluence threshold as a function 
of base temperature as determined from several series of 
data such as that in Figurc[2fb). This further emphasizes 
the softening that occurs in the insulating state with in- 
creasing base temperature. 

Summarizing the dynamics in Figure [5] and Figure [31 
(i) the conductivity rise time of ~100 ps is substantially 
longer than the excitation pulse or electron thermaliza- 
tion time, (ii) a fluence of 12.8 mJ cm~^ heats the sam- 
ple well in excess of Tc, (iii) despite (ii), the maximum 
possible conductivity is not obtained at 300 K indicat- 
ing a stiffness with respect to driving the IM transition 
and, (iv) this stiffness towards fully driving the IM tran- 
sition decreases with increasing temperature indicating 
a softening of the insulating phase. While a complete 



description of the dynamics is difficult, in the following 
we present a simple dynamic model using Bruggeman ef- 
fective medium theory that, to first order, describes the 
experimentally measured dynamics. 

Bruggeman effective medium theory (BEMT) [H, ^ 
is a mean field description to describe inhomogeneous 
media. For VO2, this corresponds to the coexistence of a 
metallic volume fraction {fm) and an insulating volume 
fraction (1-/^) which depend on temperature. Previous 
work has described the temperature-dependent conduc- 
tivity (both the finite transition temperature width and 
hysteresis) in terms of the coexistence of metallic and 
insulating phases in VO2 [l]- More recently, additional 
experimental support for this idea from time-integrated 
optical conductivit y a nd scanning probe measurements 
has been presented]!^, . 

BEMT describes the conductivity as follows: 
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where am is the conductivity in the metallic phase 
(1000 cm~^), and ai = is the conductivity in the 
insulating phase. As with previous descriptions of VO2 
using BEMT, we take the two-dimensional form (rf=2) 
of this expression ((25I ). In this simple model, there ex- 
ist disconnected metallic domains in the the insulating 
phase. Percolation of the metallic domains occurs at fm 
= 0.50, at which point the sample becomes conducting. 

We can calculate the metallic volume fraction, /,„ us- 
ing Eq. ([1]) and the experimental results presented in 
Figure[T]Jc) for temperatures above Tc (♦ for T| and • for 
T| ) , which we plot in Figure [4l The increasing temper- 
ature branch of the metallic fraction, /m(T|), increases 
from 0.52 at 330 K to 0.98 at 350 K, while the decreasing 
temperature branch, /m(T|), returns to 0.61 at 330 K, 
a consequence of the conductivity hysteresis exhibited in 
this material. In the insulating phase, we cannot use this 
approach to determine /,„ since the conductivity is below 
our detection limit. 

To describe the conductivity dynamics in the insulat- 
ing phase using BEMT we determine the temporal de- 
pendence of the volume fraction using the following ex- 
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pression: 



(2) 



With this simple model, the growth rate of fm depends on 
directly on the available nonmetallic fraction (!-/„), 
and which describes the rate at which /,„ evolves. 
It is reasonable to assume /3(T) = /3o exp(— O/fcbT) 
which describes an Arrhenius-like temperature depen- 
dence where O is an energy barrier related, in the present 
case, to the latent heat. For example, for homogeneous 



domain growth, 9 oc (T — Tc) [27|. The temperature 
dependence of I3{T) is important to consider as there is a 
bath (i.e. the substrate) to which deposited heat can es- 
cape. This heat escape, described as T = Tq exp(— t/Tgub) 
(where tau derives from the thermal mismatch between 
sample and substrate) imparts a time dependence to 
/3{T). This allows for the parameterization of (3{T) in 
terms of t. A subsequent Taylor expansion of (i{T{t)) 
about t = yields an analytical solution to Eqn. 2 
given by /™(i) = $/(! + <I>) where $(i) = /^/(l - PJQ 
where is the initial metallic volume fraction, and 
C = exp(Tj„^/35(l-ea;p(-VT-;„J). It is the term C which, 
even in the presence of superheating well above Tc, pre- 
vents the full conductivity from being obtained. The 
superscript r indicates that dimensionless factors (> 1) 
from the Taylor expansion have been incorporated into 
the effective lifetime and and rate to simplify the expres- 
sions in this phenomenological model. 

This solution describes the situation where the rate 
of increase of decreases as energy initially deposited 
in the film escapes to the substrate. This determines 
the rise time of the conductivity and maximum induced 
change which in turn depends on the initial tempera- 
ture (immediately after heating) and the initial volume 
fraction. We emphasize that this solution describes the 
rise time and must be multiplied by another exponen- 
tial exp(— t/xd) to include the subsequent conductivity 
decay (parameterized by r^). The solid black lines in 
Figure [21^ a) are fits using this description where, for two 
dimensions (j{t) = {2fm{t) — l)crm- For the fits, rj^^ ^ 
100 ps, /35 = 3.2xl0iOs-i, Td = 1 ns, and ^ 0.08 
(300 K), = 0.13 (310K), and = 0.3 (320K). The 
values of in the insulating phase arc plotted in Fig- 
ured 

We see that BEMT, appropriately extended to describe 
a dynamic metallic volume fraction, can account for the 
observed conductivity dynamics and strongly suggests 
a scenario where metallic precursors grow and coalesce 
upon photoinduced superheating. Furthermore, the re- 
sults display an initial condition sensitivity described by 
the initial volume fraction We note that our analysis 
has assumed homogeneous growth of fm from an initial 
/*„. It is possible that there is also photoinduced nucle- 
ation in which case the values of will be smaller than 
what we have estimated from our analysis. Nonetheless, 



even in the case of photoinduced nuclcation, the exper- 
imental data still reveal an initial condition sensitivity 
consistent with softening of the insulating phase and the 
BEMT describes the essence of the observed conductivity 
response. 

In summary, we studied the near-threshold behavior of 
the photoinduced phase transition in VO2. For the first 
time, we use optical pump THz-probe measurements to 
directly measure the change in conductivity of the sys- 
tem. The observed dynamics of the photoinduced phase 
transition result in the enhancement of fluctuations as 
the temperature is increased towards the transition tem- 
perature. These results may also be conducive to high- 
sensitivity optical devices, which make use of correlated 
oxides for switching, detection or optical limiters. 
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